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Abstract 

Within the ideas of pseudo-supersymmetry, we have studied a non-Hermitian Hamiltonian H- — 
a;(^^^ + |) + Q^'^ + /3^^^, where a ^ 13 and ^ is a first order differential operator, to obtain the partner 
potentials V+{x) and V-{x) which are new isotonic and isotonic nonlinear oscillators, respectively, 
as the Hermitian equivalents of the non-Hermitian partner Hamiltonians H±. We have provided 
an algebraic way to obtain the spectrum and wavefunctions of a nonlinear isotonic oscillator. The 
solutions of V-{x) which are Hermitian counterparts of Swanson Hamiltonian are obtained under 
some parameter restrictions that are found. Also, we have checked that if the intertwining operator 
satisfies rjiH- — H+rji, where t^i — p^^Ap and A is the first order differential operator, which 
factorizes Hermitian equivalents of H±. 
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1 Introduction 



A class of quantum potentials known as conditionally exactly solvable(CES) models which include one or 
more specific fixed coupling constants was introduced by Dutra [Ij . Since then, CES models have attracted 
considerable attention [31[31|31[5J[Sl[71[Slin]- Recently, a non-polynomial, one dimensional isotonic nonlinear 
oscillator potential [SI [71 |H1 IS] and its generalization to a c? dimensional model have been studied [TU] . On 
the other hand, specific symmetries of the isotonic potential have been investigated [TT]. Previously, a 
CES model which is an isotonic nonlinear oscillator potential constructed by a superpotential in [3] . The 
isotonic potential which is a harmonic oscillator with a barrier potential is given by [H [5] 

V{^) = + ^5^"' (1) 

where g > Classically, isotonic potential model has a property which is the family of periodic 

solutions with the same period |12| and quantum mechanically the energy spectrum is equidistant ^l3] . 
Recently, for the Schrodinger equation, Cariiiena et al investigated the solutions of 

where the new additional term is non-polynomial and it is the sum of two centripetal barriers in the 
complex plane [B]. Later Fellows and Smith showed that this potential is a supersymmetric partner of 
a harmonic oscillator [7]. Within position dependent mass aproach, it is shown that the source and the 
target potentials share the same energy where the isotonic potential is the reference potential [5]. At 
the same time, the fundamental framework of pseudo-Hermitian quantum mechanics has been proposed 
in [14] and the idea of non-Hermitian Hamiltonians has spread to many branches of physics such as 
electrodynamics [TS], biological physics [TS^, quantum information [17 and theoretical physics [T51 [TOl 

In this study, we discuss the isotonic oscillator with a non-polynomial term which is a Hcrmitian 
counterpart of non-Hermitian Swanson Hamiltonian [521 HSl IMl HSl HSl HZ] dHl llS]j in other words, our 
model can be a generalization of the non-Hermitian models for a nonlinear isotonic potential within the 
framework of pseudo-supersymmetry [30l IHl [32] ■ We shall also seek for the solutions of the corresponding 
system. The plan of the paper is as follows. In section 2, we introduce a brief review for factorization 
and hierarchy of Hamiltonians. Section 3 devoted to non-Hermitian partner Hamiltonians H± which are 
isospectral except for the ground state. Pseudo- supersymmetry relations between the partner Hamil- 
tonians are discussed. In section 4, the solutions of the Hermitian equivalent Hamiltonians of H± are 
introduced. 

The aim of this work devotes us to use factorization approach. As is well known, the superalgebra 
constructed by the super Hamiltonian T-Lg and the supercharges Q, is [33) 

Hs = {Q,Q^} (3) 

where Q, commute with Us and 'f is the Hermitian conjugation; 

-g^' = [Q,H.] = [g^H.] -0. (4) 

The supercharges are given by the following matrices, 

0\ ^+ /O At 
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Here, Us is diagonal that can be factorized in the manner 



'^s =i^n ^ ], n-= A^A, H+ = AA^ (6) 



H- 

where A and are first order differential operators given by a more general realization |34) 

A^a{x)^+b{x), = -a{x)^+b{x)~a'{x), ' = (7) 

ax ax ax 

where b{x) is the superpotential related to the partner potentials 

U+{x) = U-{x) + 2h'{x)a{x) - a{x)a"{x), (8) 
which are elements of the partner Hamiltonians: 

The ground states can be introduced as 

Atl!^ or A^+ = (10) 

where ■0,^ [n = 0,1,...) stands for the bound state wavefunctions for 'H±. Let Ef[ are the energy 
eigenvalues of Hi, then we suppose that the "H^ is known and has normalizable vanishing ground state 
energy in the unbroken SUSY case. In this manner, we have 

i^-{x) = {E+)-'^A^t. ^+ ^{E~)-iAi^-, E+ = E-^,, Eo=0. (11) 

2 Non-Hermitian partner Hamiltonians 

In order to describe a physical system, a non-Hermitian Hamiltonian should be 77-pseudo Hermitian 

h1_ = T]H_Tr^ (12) 

where 77 is a linear, invertible operator. If rj is not a Hermitian operator, then the Hamiltonian is 
known as weakly pseudo-Hermitian |35j . To find a Hermitian Hamiltonian h-, one can use a similarity 
transformation h- = pH-p~^ and the metric operator is obtained using 7] — p^. Consider the ladder 
operators that are first order differential operators 

^ = a{x)^+b{x), ^^a{x)^+b{x)~a'{x) (13) 

ax ax 

which will be used in a pseudo-Hermitian quadratic Hamiltonian which is known as Swanson model |22j : 

H^=LuiCk+l) + c,e+P^^\ a^f3. (14) 

where parameters w, a, /3 e M and it was shown that the eigenvalues are real and positive — 4a/3 > 
when ^ represents linear harmonic oscillator annihilation operator. Now, (jl4p can be cast in a differential 
operator form: 

= -u-^a^{x)-^ +bi{x)-^ + ci{x) (15) 
ax ax ax 
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where uj = uj — a~/3,uj>0. One gives bi{x) and C2{x) as [5S] 

bi{x) = (a - I3)a{x){2b{x) - a (x)) 

and 

ci{x) ^ {uj + a + (3)iP-{x) -{uj + 2l3)a {x)h{x) - (w - a + I3)a{x)h {x) 
+(3{a{x)a' {x) + a'^{x)) ~ 5 a {x) + ^. 

Hermitian equivalent of (|14l) which is /i_ 

obtained by using p — exp ^ / dx^^^^^ , then takes the form [151 [2H] 

d n , , d 
h- = -(D— a'^(x)— + V"_(a;) 
ax ax 



V_{x) = ( + cj + 2(a + /3) ) 5(x)(&(x) - a {x)) - (w + a + ;3)a(a;)& (x) 



where 



+ -^—a{x)a (x) + - I + 2{a + /3)\a (x) + . 

In order to factorize , it is written as a product of two operators 

Let us give an ansatze for the operators we use in above 

A — V^a(x)— h b(x), A^ = —V^a(x)— h b(x) — V^aix). 

ax ax 

where b(x) and a{x) are real and they can be chosen to generate V±{x). Thus, V-{x) becomes 

V^ix) = P{x) - V^^{aix)b{x)). 

ax 

As is well known, /i_ is isospectral to its partner Hamiltonian h-^ that is 

h+ = AA^ 

and it can be expressed as 

h, = -^±a^x)^^+V,ix). 
Now, the partner potential of V- (x) is 

V+(x) — o{x) + \fGj I —\fG)a{x)a(x)" + a{x)b' (x) — b{x)a [x) 



Exploiting the fact that the non Hermitian partner Hamiltonians H± and their Hermitian counterparts 
h± are related by a similarity transformation 

H± = p-'h±p, (27) 

then we can get ff+ as 

= -Co^a\x)^ + h,{x)j- + I - + V+{x). (28) 

ax ax ax I 4a;a^ 

Finally, we can give a diagramatic representation for a chain of Hamiltonians we used as follows 

^ ^ h+ ^ H+ (29) 
such that we can arrive at partner Hamiltonian using this chain. 

2.1 Factorization of the Hermitian counterparts of Swanson Hamiltonian 

Some choices of a{x) and b{x) may lead to an effective potential model. In (331 [5S], it was shown that 
when the commutator is constant, i.e [^,^^] = 1, then the relation between a{x) and b{x) can be given as 
[281 



The more general case corresponds to [^,^^] 7^ 1. Then, let us use the ansatz for a{x) and h{x) 

a{x) = x^ (31) 



in (PO)) . then V^{x) becomes 

V-{x) = ^ + 2(03 + 204)2;^ + (-2ai +a2)(c+ 1) + 05 (33) 
x^ 

(2ai + oic + 2aid — 3a2d)x^ + 2ai(l + d) — a2d 
(.t2 + d)2 

where c, d are some real constants. Here a^, i = 1, 2. ..5 can be recalled from (PU)) . i.e: 

ai = ^"^^^ +lj + 2(q + /3), a2=£l' + a + /3, 03 = ^^i^ (34) 
cj 2 

04 = -A ^ J^' +2(Q + /?)j, 05 = (35) 

With an eye on what is to follow, let us first match ([55)1 with ([3^ to factorize so another ansatz is 
introduced for h{x) which is, 

bix) — — — Qx ^ — T-—,- (36) 

^ ' X x^ + d ^ ' 
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At that case, V- (x) is written as 



fi^ ^ , -(2g + v^)a:4 + {-'id^ + A + 2m - 2dg)x'' + 



Hence, we can give the couphng constants of V-{x) in (j37p in terms of uj,a,f3 so that we can compare 
([37l) and ((33|). At the same time, its partner V+(x) becomes 

V+[x) = -^ + (g +qVu}-2u})x -/x(2g+3Vw)+A — ^-15 . 

(38) 

After some straightforward algebra, partner potentials happen to be 

V- {x)^^ + gig + 'SV^)x^ ^{h + X){Vq + 2g) (39) 
x^ 

{2gd + 2^ + \- d^)x^ + d(2/i + d^/^ + 2gd) 

V+{x) ^ t-^ + g{g + - 2u))x^ - (a* + A)(3\/^ + 2(?) 
(2g(i + 2^ + A + 5dv^)a;2 + d(2/i + + 2gd) 

We seek for a solvable V+(x) which implies some parameter restrictions such that the non-polynomial 
part in (38) vanishes within these conditions below 

A = -2dV^, M = -^(2g + 3V^). (41) 
Herefrom, if we plug (|4T|) into V+{x), it turns into 

V^ = l^ + ig^ + gV^- 2C)x^ + 2d^g+ I^) (^^ + ^) • (42) 
The same procedure is followed for V- (x) , it follows that 

y--i + 9i, + 3V^).^ +2d(g+ I^) + ^) + ^^d'ff^- (43) 



Comparing ([55)1 and as we stated before, it may give us a chance to write /i. A, g in terms of w, a, {3. 
Then, 

ai = = —{2g + i^f (44) 
2(03 + 2a4) = £-(g + 3v^) (45) 
(a2-2ai)(c+l) + a5 - 2d(g + I^)(g + ^) (46) 
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and 



c{~3a2d + 2ai + aiC + 2aid) = 12LU(f (47) 
2ai{l + d) - a2d = Auid^ . (48) 

We can solve c using last two relations above: 

ci2 = ( 2a2d - iidd^ ± JUojd^ - 2a2dV + iSaiCudA . (49) 
2ai \ ) 

Using (gH), (|45]) and (|46]), q can be found as 

Q = ^(±\/4X - 27w - 4\/^) (50) 

where 

^ = - 8(03 + 2a4) + 7^ ((a2 - 2ai)(c + 1) + ag) . (51) 

2a 

There are some restrictions on parameters such as > 0, so we take q — ^{\/AX — 27(D — 4^/5). In- 
side the square root in ((50)) must be positive, g is positive, then AX > ASoj. Because 02 — 2ai — 
_ ui +3{a+i3)uj+2{a-0) ^ then wc need c < 0, |c| > 1 and we use the negative sign for c which is 
c = ^{2a2d- Aujd^ - ^J{Au}d^ - 2a2dp + ASaiuid^). When ((SO]), (011) are satisfied, then Hermitian 
counterpart of iJ_ which is /i_ given by ([T9| is factorized. 



2.2 Pseudosupersymmetry 

Here, is diagonalizable with a discrete spectrum which admits a set of complete biorthonormal 
eigenvectors {|*n), |*n)}: 

i/_|*„> = ifl I $„> (52) 

^|<l>„)(*„| = 5]|vl'„)(<l>„| = l. (53) 

n n 

On the other hand, the intertwining relations between h± can be written as 

h-A^ = A^h+ (54) 

and 

h+A = Ah^. (55) 

And, given in (jl4p is a non-Hermitian and diagonalizable Hamiltonian with real or complex conjugate 
eigenvalues. If there exists an operator rji such that 

mH- = H+r^i (56) 

where is the partner Hamiltonian of then, the intertwining operator rji is given by r/i = p^^Ap 
where A and its adjoint satisfy ([51]) and ([55]) . This relation 771 = p^^Ap was proven before [SIJ- The 
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intertwining operator that links a non-Hermitian Hamiltonian to the adjoint of its pseudo-supersymmetric 
partner Hamiltonian was also studied |32) . If one uses the similarity transformation 

H+ = p-'h+p (57) 

then, if we multiply (j57p by p~^A from right, we have: 

H+p-^A^ p'^h+A (58) 

and if we use ([55]) in ([55]) . we get 

p-'^Ah^ = H+p-^A (59) 

p-^App^^h^p = H+p-^Ap^ (60) 

p-^ApH^ = H+p-^Ap. (61) 

Then, it is seen that the intertwining operator 771 is given by 771 = p^^Ap. Following the same way we 
can obtain 772^?+ = H^ri2 where 772 = p^^A^ p. Here, 771 and its pseudo-adjoint 772 can be shown as 

vl = v'^vlv = P^^-^P = m (62) 

which leads to construct pseudo- super algebra of non-Hermitian supersymmetry. Thus, the operators 
Q, become 

The pseudo- super Hamiltonian can be introduced as 

^-{%- h) (") 

where pseudo- super charges satisfy 

H = [Q,Q«]+ = 0. (65) 
Now, let us find rfi for our problem, it can be written as 

771 = V^x'^ — + _ ^ , -{{a-f3-gV^)x^ + {{a^l3){d-c-l)-{^ + 2gdV^))x^+d{pV^-a + (3)) 

dx y'ujx[d + x'^) 2 

(66) 

which connects two partner Hamiltonians H±. We can arrive at using 771 given above. 

3 Solutions 

In this section, we aim to obtain the solutions of the partner Hamiltonians with (j42|) and (|43|) . Let 
eigenfunctions of h± and h± be ^p^{x) and if>^{x). Then, h± can be transformed into h±, 

h± = eh±9-^ = -ua^{x)^ ~ Cja{x)a'ix)4- + V±ix) + -a'^{x) + -a{x)a"{x) (67) 

dx'' dx 4 2 
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where 6 = a{x). One may use the change of independent variable in (|67p as, 

dy 



/uja{y) 

and for our problem it equals to z = — i^. Then we arrive at 



where we can write partner potentials V± — V± + ja'^{x) + ^a{x)a" [x) in terms of z as 



and 



(68) 



h± = ~j^+ {V^{x) + ja'-\x) + '^a{x)a"{x))^^^ (69) 



y+(z) = /iw2;+ — + 2d[Q + ^^){q + (70) 



T~r / X 9-9 g^ + Sv^e + 2a; 1 „ 7^^J^/ , _ , , _ , 2da;z^ — 1 ,„_, . 

w 2 2 (dwz^ + \y 

We know that ([70]) and ((7T|) are isospectral. In [36], an eigenvalue equation which is given by 

- — + — + 5z2 j x„(z) = e„X„(z) (72) 

has exact eigenvalues and eigenfunctions as 

e„ = 2(5(2n + 7), n = 0,l,2,... (73) 

where 7 = 1 + iyTT^J, 5 = \fB and 

X„(z)=iVz^-i/2e-*^'/2 iFi(-n;7;<5z2), 0,1,2,... (74) 

This model was studied by Goldman and Krivchenkov |3 7] who showed that the energy spectrum of 
this potential is an infinite set of equidistant energy levels. Then, if we use uj — ^^(2p + 3-^/5) and 
7 = + \ for the potential V+(z), exact eigenvalues and eigenfunctions of V+(z) are given by 

E+ ^ 2u}{2n + ^ + 5) (75) 



/UJ 

and 

^tiz) = C„z^-^e-^"^' iFi(-n;7;a;z2) (76) 
where C„ is the normalization constant which is given by |36j 
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and the solutions (^^(2) are square integrable on (0, 00) [35] ■ The relation between the confluent hyper- 
geometric function and Laguerre polynomials may be used later: 

,F,{-n;b+l;y) = ——-Lliy). (78) 

(b+l)n 

Here the Pochhammer symbol (s)„ = ^p^"^"^ is used for both ([77]) and ([78]). Now, in order to find the 

eigenfunctions of /i_ we can follow an algebraic way such that we can start with a mapping between h± 
and h± introduced before: 



Also, this relation can be written as 



h± = 0h±e-^. (79) 



h+ = OAA^e-^ (80) 
= AA^ 



and following the same way one obtains 
where A and A'' can be introduced as 



= A'' A (81) 



A = V^a{x)-^+b{x)-—a'{x) (82) 
ax 2 

A^ = -^a{x)^+b{x)-—a'{x). (83) 
dx 2 

Finally we can use A^ to find <y3^(z) which is the solution of V-{z): 

(p-{z) = = (^-V^a(x)^ + b{x) - ^a'{x)^ C'^z^-^/^e-'^^'^' Ll-\uz^) (84) 

f d pi 2dLdZ \ r „ 1 /o l,-,,2 ,,1,9, , , 

= [--^+u:z + ^-+ A C„z^-i/^e'^- Ll'\(,z^) (85) 

\ dz ^Lo z 1 + dwz^ / 

where we use /i and A from (|41l) in (|84l) . Using a relation [TU] 

|l^(0 = -L^ll(t) (86) 

and the identity of Laguerre polynomials 

^,{t)^Li_,+Li-\t) (87) 

helps us to re- write <p^(^) as |10j 
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On the other hand, one can show that the solutions ipn{z) can be normaHzed. For example, ip^ can be 
given as 

ib+iz)^Afz'>+^^^e-i^'" iFi{-n;r,^z^) (89) 
where Af is the normahzation constant. To normahze t/j^, we may introduce the integral [3 8) 



oo 



(90) 



for u > 0, m, n = 0, 1, 2, ... and dmn = if to 7^ n, Smn ~ 1 ii m = n. Now we use 
then we have 

hh (7)^(7).fc!^! Jo ^ ^ 

The integral representation of the gamma function is 

T{x) = / e'H^-^dt, x>0, (93) 
Jo 

and using a variable change, we obtain 



n / m. 



J _ 1 V- (-TO);(7 + fc + l)i (-n)fcr(7 + fc + l) 

^™ - 2^7+1 1^1. (7).^! J (7).^! ■ ^ ^ 

The term in the bracket in above relation corresponds to hypergeometric function 2^1, then we re- write 

Jmn ^S 

■'••"^-^t .^(-..;7 + * + i;7;i) '-'-;!7:,^"" ^ M 

fe=o ^''^ ■ 

Using Chu-Vandermonde identity which is given by |36| 

r(c)r(c-a-b) 

2i^i a; 6; c; 1) = — — (96 

1 (c — a)i (c — b) 

T - ^(^) A (-n)fer(fc + 7 + l) r(m-fc-l) 

^,nn 2c2;7+lA. r(7 + TO)r(-fc-l)- ' 

Using the identity [35] 

r (-i)"/c! n < 77 < A" 



in (|95)) . we get 



and (|97p . we obtain 
for k = n = m. Finally, normalization constant becomes 



_ (n + 7)(n + l)!r(7) 
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4 Conclusion 



To conclude, we have derived a new class of isotonic nonlinear oscillator using the concepts of pseudo- 
supersymmetry. We have studied a generalized non-Hermitian Hamiltonian and its partner that can 
generate solvable isotonic and nonlinear oscillator potentials by using an appropriate mapping, i.e. a 
similarity transformation. Hermitian partner Hamiltonians which include effective solvable new isotonic 
and its nonlinear partner potential are factorized and it is seen that the nonlinear oscillator is solvable 
under parameter restrictions given by (1411) . We have also checked that the intertwining operator rji 
that connects non-Hermitian partner Hamiltonians H± can be given in terms of a transformation of the 
factorization operator ^ as 771 = p~^Ap. At the end, we give solutions of isotonic nonlinear oscillator 
after a coordinate and eigenfunction transformations. We have seen that the eigenvalues are real and 
positive which agrees with the results of We introduce formulae for all the discrete eigenvalues and 
normalized eigenf unctions. 
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